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I. INTRODUCTION 



One ultimate goal of contemporary strong interac- 
tion physics is to find a comprehension of the physical 
properties of hadrons by means of the underlying the- 
ory of Quantum Chromodynamics (QCD). Among sev- 
eral methods which provide a link between QCD (quark 
and gluon) degrees of freedom and the hadronic spec- 
trum are the QCD sum rules which have to be considered 
as important nonperturbative approach in understand- 
ing the physical observables of hadrons. The sum rule 
method, first developed for the vacuum [lj, has later been 
extended to finite density |2|. y. finite temperature 
|M and mixed finite density and finite temperature 
7J. Within the QCD sum rule approach, and more gen- 
erally in hadron physics, pions and nucleons have to be 
considered as important degrees of freedom because the 
pion is the lightest (Goldstone) meson, while the nucleon 
is the lightest baryon. In-medium QCD sum rules pro- 
vide a direct way to relate changes of hadronic properties 
to changes of the various condensates, i.e. nucleon and 
pion expectation values of quark and gluon fields. There- 
fore, expectation values of a local operator O taken be- 
tween these states, (7Tp hys |0|7rph ys ) and (A^ phys |0| N phys ), 
need to be known. However, the predictive power of 
the QCD sum rule method in matter meets uncertain- 
ties when evaluating condensates, especially higher mass 
dimension condensates inside the nucleon. Accordingly, 
the exploration of nucleon matrix elements is presently 
an active field of hadron physics, cf. H,0|- 

If the operator O consists of hadronic fields, then in 
principle one needs an effective hadronic theory which 
decribes the interaction between pions and nucleons, re- 
spectively, and the hadrons from which the operator O is 
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made of for evaluating these matrix elements. However, if 
one is concerned with pion matrix elements then the use 
of soft pion theorems [T(J, HH Il2> LLj 13 gives in general 
good estimates for such expressions, which are related 
to several so called low-energy theorems like Goldberger- 
Treiman relation , Adler-Weisberger sum rule | 16| or 
Cabibbo-Radicati sum rule 17]. These soft pion theo- 
rems as algebraic tools are based on the hypothesis of par- 
tially conserved axial vector current (PCAC) 0,0,|2(j 
and postulated current algebra commutation relations 
|l2l Il3| , and allow in general to trace the pion matrix 
elements of operators made of effective hadronic fields 
back to vacuum matrix elements. A feature of the soft 
pion theorems is that they can also be deduced within 
quark degrees of freedom. Accordingly, pion matrix ele- 
ments of quark field operators have also been evaluated 
by means of the soft pion theorem (if we speak about 
the soft pion theorem then we mean the special theorem 
considered in the Appendix ^ which is the relevant one 
in our context) expressing the pion field and axial vec- 
tor current, respectively, by interpolating fields made of 
quark degrees of freedom 0, [2^] . 

After discovering the powerful method of current al- 
gebra for mesons several attempts have been made to 
investigate the possibilities for extending this algebra to 
the case of baryons. Especially, the analog hypothesis of 
a partially conserved baryon current (PCBC) and the re- 
lated (and postulated) bar yon current al gebr a has been 
investigated long time ago |2ll2l I2I El I27L Esf . These 
attempts focussed on the construction of baryon currents 
by products of nucleon fields. Furthermore, in |29| this 
procedure has been studied by considering baryon cur- 
rents made of quark degrees of freedom where several 
relations between form factors, e.g. baryon-meson ver- 
tex form factors, have been obtained. However, it turned 
out that, while the PCAC directly leads to the mentioned 
soft pion theorems for evaluating pion matrix elements, 
the PCBC does not provide a comprehensive algebraic 
theorem for evaluating nucleon matrix elements. There- 
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fore, up to now for evaluating nucleon matrix elements 
of an operator O consisting of quark fields more involved 
tools are needed like chiral quark model 111, lattice evalu- 
ations , or Nambu- Jona-Lasinio model[3(j ■ From this 
point of view it seems very tempting to look for an al- 
gebraic approach for evaluating nucleon matrix elements 
in analogy to the soft pion theorem. Here, by using di- 
rectly the nucleon field instead the nucleon current, we 
propose such an algebraic approach for evaluating ma- 
trix elements of quark operators taken between a bare 
nucleon, i.e. the valence quark contribution. 

To clarify what the terminology "bare nucleon" means 
we recall the basic QCD structure of nucleons. From deep 
inelastic lepton-nucleon scattering (DIS) experiments we 
know that nucleons are composite color-singlet systems 
made of partons. In the language of QCD these are three 
valence quarks with a current quark mass, accompanied 
by virtual sea quarks and giuons. Accordingly, the phys- 
ical nucleon state |iVph ys ) is a highly complicated object 
consisting of many configurations in the Fock space. For 
instance, in the case of the proton, the Fock expansion be- 
gins with the color-singlet state \uud) consisting of three 
valence quarks which is the so called bare proton state, 
and continues with \uudg), \uudqq) and further sea quark 
and gluon states that span the degrees of freedom of the 
proton in QCD. 

In the low energy region, many properties of the nu- 
cleon can rather successfully be described by approxi- 
mating the virtual sea quarks and giuons by a cloud of 
mesons, especially pions, surrounding the bare valence 
quark core. Accordingly, in the pion cloud model, which 
resembles the Tamm-Dancoff method [U IH |H |H 
[3|| the physical nucleon is viewed as a bare nucleon, 
which accounts for the three valence quarks, accompa- 
nied by the pion cloud which accounts for the virtual sea 
quarks and giuons. Then the Fock representation for the 
physical nucleon reads plEs! I37tl38l Hoi] 



|iV phys ) = Z'J 2 (\N) + <\> X |iV7r} + 2 \NtT7t) + ...) , (1) 



where the Fock state | N) represents a bare nucleon state, 
|iV7r) and | A^7T7t) represent a bare nucleon with one pion 
and two pions, respectively, and the dots stand for all 
of the Fock states consisting of one bare nucleon with 
more than one pion or heavier mesons. The probability 
amplitudes <j) n to find the nucleon in the state |iVrar) 
can be evaluated by using a Hamiltonian which describes 
the pion-nucleon interaction [H El HI IM E3 • Then 
the bare nucleon probability can also be determined and 
turns out to be Zn — 0.9 |32ll3S| . Since the deviation of 
Zjy from 1 comes from pion-nucleon interaction one has 
to put Zn = 1 if the pion cloud is not taken into account. 
By using the Fock expansion the expectation value 
of an observable O taken between the physical nucleon 



states is given by [H |U , 

(N phys \6\N phys ) = Z N (^(N\6\N} + 

<t>\ {Nt:\6\Ntt) + c/4 (Nirir\6\Nirir) + ...j . (2) 

The first term on the right side of (J5J, i.e. the contribu- 
tion of the bare nucleon without pions, plays an impor- 
tant role for two reasons. First, the bare nucleon is ex- 
pected to give the main contribution in many cases |42| . 
And second, for the leading chiral correction one needs 
only the contributions of the lowest-momentum pions in 
the cloud allowing an application of the soft pion theo- 
rem (see Appendix^}, which then reduces the pion cloud 
terms in also to bare nucleon matrix elements [43Ll44| . 
Accordingly, in this paper we focus on bare nucleon ma- 
trix elements and propose an algebraic method for eval- 
uating them. This approach seems capable to estimate 
nucleon matrix elements of quark operators in a straight- 
forward way. We also note that within the algebraic ap- 
proach new parameters are not necessary since the bare 
nucleon matrix elements are traced back to vacuum ma- 
trix elements, like in the soft pion theorem. We apply 
the method on two-quark, four-quark and, finally, on six- 
quark operators inside the nucleon which so far have not 
been evaluated. 

The paper is organized as follows. In section^ we de- 
rive an algebraic formula for evaluating matrix elements 
taken between the state of a bare nucleon. In section UTTI 
a valence quark field operator with the quantum numbers 
of a bare nucleon is introduced. A few tests of the nucleon 
formula on well known bare nucleon matrix elements of 
two-quark operators are given in section HV Al (currents) 
and IIV Bl (chiral condensate) . In section IIV CI we ex- 
plore the valence quark contribution of four-quark con- 
densates within the algebraic method developed and as- 
sert an interesting agreement with the results of ground- 
state saturation approximation when taking properly the 
valence quark contribution. We also compare our findings 
for the valence quark contribution of four-quark conden- 
sates with recently obtained results within a chiral quark 
model. In section we evaluate six-quark condensates 
inside the bare nucleon. A summary of the results and 
an outlook can be found in section IVT1 In Appendix lAl a 
derivation of the soft pion theorem is given which shows 
the similarity of it with our algebraic approach. Details 
of some evaluations are relegated to the Appendix FBI 

II. NUCLEON FORMULA 

Let 0(x) be a local operator which may depend on 
space and time, x — (r,t). We are interested in matrix 
elements taken between two bare nucleon states \N(k, a)) 
with four-momentum k and spin a (i.e. \N) is either a 
bare proton \p) or a bare neutron \n) state, which are con- 
sidered as QCD eigenstates). To derive a formula for such 
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matrix elements between bare nucleons with finite nu- 
cleon masses and momenta we first apply the Lehmann- 
Symanzik-Zimmermann (LSZ) reduction pH lifij on one 
nucleon state, 

(7V(fc 2 ,a 2 )|0(x)|iV(fc 1 , ( 7 1 )} 



iZ, 



-1/2 



d 4 Xl {N(k 2 ,a 2 )\ T w (6(x) i a N (xi)^j |0> 



x f «7m 9^ + M N ) u^(fci,(7i) 

V / a (3 



-ik\x\ 



(3) 



where the greek letters a,/? are Dirac indices. The 
normalization of the nonperturbative QCD vacuum is 
(0|0) = 1, and the normalization for the nucleon 
state reads (JV(fc 2) g 2 )|iV (fci, oi)) = 2£ fel (27r) 3 <5( 3 )(ki - 
1*2 ) S aia2 , where E^ — ^/k^ + M'fj. Throughout the pa- 
per we take the sum convention: If two Dirac (or later 
color) indices are equal or not given explicitly, then a sum 
over them is implied. The four-momenta are on-shell, 
kf = k 2 = ^n-i f° r noninteracting nucleons the bare nu- 
cleon mass equals the physical nucleon mass, M/v = 938 

MeV. The field &n(xi) is the interacting (adjoint) nu- 
cleon field operator, i.e. off-shell. The equal-time anti- 
commutator for the interacting nucleon field operator is 
the same as for the free fields and reads 



= <5^(n -r 2 ) S a ^ 



(4) 



For the wave function rcnormalization constant we have 



Z' 



-1/2 



< f. The free nucleon spinor satis- 
0, 



A(x 1 )B(x 2 )Q(h - t 2 ) + B{x 2 )A(x 1 )Q{t 2 - tt), has been 
taken in Eq. © . 

We approximate Eq. © by introducing a noninteract- 
ing nucleon field operator given by 



d 3 k 1 

(2tt)3 2E k 



^2 (a N (k,a) u%{k,a) e~ 

cr=l ^ 

+P N (k,a) v%(k,a) e ik *) ,(5) 



with the corresponding anticommutator relations in mo- 
mentum space 



a N (ki,(Ti), a' N (k 2) <j 2 ) 
= 2E kl (2tt) 3 £( 3 )(k! -k 2 )<5 CTl 



b N (ki,ai),b* N {k 2 ,a 2 ) 



+ 
(6) 



Accordingly, \N(k,a)} — a' N (k, <t)|0). For the noninter- 
acting nucleon field operator Z^ = f , and the equation 
of motion follows from (J5J, (ij^ d% — Mj^j ^n{x) = 0, 
and for the adjoint noninteracting nucleon field operator 
it reads ^n(x) {id^^ + M/v J = 0, respectively. Then 
one arrives at 



(N(k 2 ,(72)\0{x)\N(k 1 ,a l ))= I d 4 Xl e-^ S(t-h) 



G(x),H> N ( Xl ) 



1°) (7o) Q/3 vr N {ki,ai) 



(7) 



fies (j^k^ — Mjv)ujv(A;, c) = 0, with normalization x (N(k 2 ,o 2 ) 
UN(k,<r 2 )uN{k,(Ji) = 2Mjv S ai(T2 . The operator O is, 
for physical reasons, assumed to consist of an even num- 
ber of fermionic fields, i.e. a bosonic operator, accord- 
ing to which the Wick time-ordering, Tw A{xi)B(x 2 ) = Applying this procedure on the left nucleon state yields 

I 



(N(k 2 ,a 2 )\0(x)\N(ki,ai)) = / d 4 Xl / d 4 x 2 e 



X < 2 (fc 2 ,Cr 2 )(7o) ;32Q2 (0| 



-ik\X\ ik^X2 



S(t-h) 5(t-t 2 ) 



d(x) 7 ^( Xl ) 



1°) (7o) ai Pl u^ik^ax) , 



(8) 



r 



which is symmetric under the replacement 

(0|[*,[d,t]_]+|0> -> (O|[[*,0]_,t]+|O>. Eq. © 
is the central point of our investigation and we call 
it nucleon formula. This formula resembles the soft 
pion theorem given in Appendix ^ ft is worth to 
underline that, due to the S- functions in JSJl, only the 
equal-time commutator and anticommutator occur. The 
anticommutator comes into due to the fact that the 
commutator in J7J) between the operator O (consisting 
of an even number of fermionic operators) and the 



(adjoint) fermionic field operator '5 /v yields an operator 
consisting of an odd number of fermionic field operators. 
Therefore, when applying LSZ (cf. Eq. ©) on the other 
nucleon state a Dirac time-ordering, TdA(xi)B(x 2 ) — 
A(xi)B(x 2 )@(t 1 -t 2 )-B(x 2 )A(x 1 )Q(t 2 -t 1 ), is needed. 

The nucleon formula is valid for a noninteracting nu- 
cleon with finite mass Mat and finite three momentum k, 
and in this respect it goes beyond the soft pion theorem, 
which is valid for pions with vanishing four-momentum 
only. In the next section we supplement the nucleon for- 
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mula with a nucleon field operator expressed by quark 
fields, which allows then the algebraic evaluation of bare 
nucleon matrix elements of quark operators. 

We note a remarkable advantage of the algebraic ap- 
proach. The operator O is a composite operator, i.e. 
a product of field operators taken at the same space- 
time point. As it stands, such a composite operator 
needs to be renormalized. Therefore, a renormalization 
& cn = 6- (6)0 (we abbreviate (6) = (0|O|0», which 
applies for products of noninteracting field operators, has 
to be implemented 46]. However, the term (O)o is a c- 
number and, according to Eq. ©, does not contribute 

because of [(0)o,^n]- = 0. Another kind of renor- 
malization for products of interacting field operators is 
also based on subtracting of c-numbers (so called renor- 
malization constants) which vanish when applying the 
commutator in the nucleon formula. Therefore, one may 
consider the composite operator 0(x) in Eq. ijSJ) as a 
renormalized operator. This feature is also known within 
PCAC and PCBC algebra, and in particular within the 
soft pion theorem. 



III. CHOICE OF NUCLEON FIELD OPERATOR 

The nucleon formula JHJ can directly be applied on a 
local operator O which consists of bare nucleonic degrees 

of freedom, e.g. O — * n & n , ^ n 7/j * n , etc. However, 
we are interested in operators basing on quark degrees 
of freedom, e.g. O = qq,qjaQ, etc. To make the rela- 
tion JHJ applicable for such cases one needs to decompose 
the bare nucleon field operator ^ # into the three valence 
quarks, yielding a composite operator ippj to be specified 
by now [42|. Although there has been considerable suc- 
cess in understanding the properties of nuclcons on the 
basis of their quark substructure as derived within QCD, 
a rigorous use of QCD for the nucleons is not yet in reach. 
Therefore, in order to gain a nucleon field operator which 
shows up the main features (quantum numbers) of the 
bare nucleon a more phenomenological approach on the 
basis of the quark-diquark picture of baryons 0] is used. 

To be specific we consider the proton. The bare pro- 
ton state \uud) is defined by the SU(2) flavor, SU(2) 
spin wavefunction of three valence quarks. In the quark- 
diquark model of the bare proton two of these valence 
quarks are regarded as a composite colored particle (di- 
quark) which obeys the Bose statistic and which has a 
mass of the corresponding meson (e.g. for QCD with 
N c = 2 Pauli-Gursey symmetry |49j), i.e. a mass which is 
significantly larger than the current quark mass. Within 
quark degrees of freedom the general expression for such 
a diquark can be written as [23, E3 

^ g2 (x) = qt T (x)Crq*(x). (9) 

Here, qf and q!? are quark field operators of flavor u or d 
with color index a and b, respectively. Throughout the 



paper all quark field operators are solutions of the full 
Dirac equation, (ij^D^ — m q )q(x) = with = — 

ig s i^A Q /2, where i° are the gluon fields and Tr(A Q A b ) = 
28 ab (a,b = 1, ...,8 are Cell-Mann indices, which should 
not be confused with the color indices (in roman style) 
a, b,c (later also i...,n) of quark fields). The equal-time 
anticommutator for these interacting quark fields is the 
same as for free quark fields and reads 

[C(r 1 ,t),^ t (r 2 ,<)] + = (n -r 2 ) 5 a p S ah . (10) 

The charge conjugation matrix is C = 17072, arid T = 
{ 1, 75, 7 M , 7^75, a^y } is an element of the Clifford alge- 
bra. C changes the parity of V, e.g. C75 has positive 
parity. The diquark JSJ, considered as a composite op- 
erator made of quark fields, does generally not commute 
with quark field operators. On the other side, if the di- 
quark is regarded as an effective boson, it commutes with 
the fermionic quark fields. This feature of the diquark, 
considered as a bosonic quasiparticle, can be retained 
on quark level when neglecting the quantum corrections 
for the quark fields which are participants of the diquark. 
Accordingly, the diquark is separated into a classical part 
and a quantum correction 

= 9i» C r q h 2 (x) + A^ q2 (x) . (11) 

The classical Dirac spinors qf, q\ are solutions of the 
full Dirac equation (ij^D^ — m q )q(x) = 0. The clas- 
sical part in i|ll|). = qf T CT q\, commutes with 
quark field operators. To specify the diquark relevant 
for a proton we note that there are only two struc- 
tures, r = 75 and r = 7570, which have positive par- 
ity and vanishing total spin, J p = + 50]. This is 
in line with |5l| . where it was found that the proton 
has indeed a large overlap with the interpolating field 

rj p = e abc ^M aT C75 d 1 ^ w c , where e abc is the total anti- 
symmetric tensor. We also remark that in lattice calcu- 
lations the field fj p is usually used [52I l53l] since this in- 
terpolating field has an appropriate nonrelativistic limit. 
In addition, the field r) p is also a part of the so called 
Ioffe interpolating field, which for the proton is given by 

ry Ioffc = 2e abc ((u aT Cd b ) 75 ?i c - (u aT Cj 5 d^)^ 0. 

The Ioffe interpolating field is usually used in QCD sum 
rule evaluations; for a more detailed motivation see also 
|55|, These properties in mind we take guide 
for constructing a proton field operator and obtain a 
scmiclassical interpolating proton field by neglecting the 
quantum correction of the diquark. Further, we assume 
that any quark of the nucleon can either be a participant 
of the diquark or can be located outside the diquark. In 
this line only two different structures for a semiclassical 
interpolating proton field may occur and, according to 
this, the general semiclassical field operator for a bare 
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proton is a linear combination of both of them: 



a b c 



A p (u aT (x) C 75 d h (x)) ii ca (x) 



+B P (w aT (x) C75 u h (x)) d ca (x) 



(12) 



The colorless operator (I12[) leads to the quantum num- 
bers of a proton (charge, parity, spin, isospin). In the 
following we evaluate proton matrix elements on quark 
level by means of the nucleon formula (jSJ) where the field 
operator *& p (x) is replaced by ipp{ x ) given in Eq. 112(1. 

Before going further we have to comment on the nor- 
malization of the field operator (|12J1 . i-e. on the deter- 
mination of the coefficients A p and B p which, in general, 
are complex quantities. With the relativistic normali- 
zation for the nucleon state and taking into account that 
there are two u quarks inside the proton we demand j5|J 
(p(k,a 1 )\u\ x ^O)u i ol {0)\p(k,a 2 )) = ±E k 8 airJ . 2 . Evaluating 
this term using the nucleon formula JSJ with 112|l we find 
the normalization, at x — 0, to be 



2 £ a kc £ a k' c ' fi CC ' y^ 1 ' ^ 



C l5 d 



,a T 



C75 d h = 2 . (13) 



From <p(Mi)l4 t (0)4(0)|p(fc>°*)> = 2E k <W 2 we de- 



duce 



\B, 



|2 ^^bc ^a b c' 



C lb u h = 1 .(14) 



Formula JSJ in combination with the field operator 112|) 
and the normalizations l|13[) and (|14|) summarizes our 
propositions made for obtaining bare proton matrix ele- 
ments for quark operators. 

For neutron matrix elements in Eq. JSJ) we have to 
insert the semiclassical field operator for the bare neu- 
tron, ip n (x), which is achieved from l|12|) by interchang- 
ing u *-¥ d, u <-> d and by the replacements A p — > A n , 
B p — > B n . The corresponding normalizations for the bare 
neutron field operator, i.e. the determination of A n and 
B n , are obtained from Eqs. I|13[) and (|14|l by interchang- 
ing the up and down quarks, and A p — > A n , B p — ► B n . 



IV. TESTING THE NUCLEON FORMULA 

In the following we will test the outlined formula, 
Eq. (JSJ with field operator Eq. (|12fl . and compare with 
known bare nucleon matrix elements. Throughout the 
paper we evaluate matrix elements of a composite oper- 
ator 0{x) at x = and therefore omit the argument x in 
matrix elements. 



A. Electromagnetic and axial vector current 

The electromagnetic current for the noninteracting 
pointlike neutron on hadronic level is zero, due to the 



vanishing electric charge of the neutron. For the non- 
interacting pointlike proton it is given by J° m (a;) = 

e p ^> p(x)t p'fy P (x) , where the electric charge of proton 
equals the elementary electric charge, e p — e. Now there 
are two possibilities to evaluate such a matrix clement 
on hadronic level: either by means of the algebraic ap- 
proach Eq. (JSJ) and the anticommutator relation JI}, or 
the usual way by means of the field operator Eq. (J5J 
and the anticommutator relations © . In both cases it is 
straightforward to show that {p(k2, a 2 )\ J^ m \p(kx, ax)) = 
e P u P (k 2 , a 2 )^f fl u p (kx, (Ti) on effective hadronic level. As 
a first test of the nucleon formula we verify this relation 
on quark level where the electromagnetic current is given 

by J* m (x) — I e u(x) 7^ ii(x) — ^ed(x)^ p d(x). Indeed, 
by using Eqs. (|B5|I and i|B6l) from Appendix [B] for the 
bare proton we get on quark level 



(p(k 2 , a 2 )\J? l nl \p(k 1 , ai)) = e p u p (k 2 , a 2 ) 7 M u p (ki, ax) 



(15) 



Similar, for the bare neutron, with Eqs. i|B7(l and 
l(B8(l from Appendix [BJ we obtain on quark level 
(n(k2, a 2 )\ J® m |n(fci, ax)) — 0. Both of these findings are 
in agreement with the results on effective hadronic level. 

Now we look at the axial vector current which on 
hadronic level for a noninteracting pointlike nucleon is 

given by A°(x) = g\ ^> n{x)~i^^ N (x) 59], where g\ 
is the axial charge of a bare nucleon. For the moment 
being in this paragraph up to Eq. (|16fl a = 1,2,3 are 
isospin indices, and ^n, \N) and un are isoduplets. The 
isospin matrices r a coincide with Pauli's spin matrices 
with normalization Tr(r a r b ) = 2S ab . 

Similarly to the case of electromagnetic current, there 
are two possibilities to evaluate this matrix element 
on hadronic level: by means of the algebraic approach 
Eq. JSJ) and the anticommutator relation or di- 
rectly by means of the field operator Eq. (J5J and the 
anticommutator relations ©. In both cases one ob- 
tains on effective hadronic level the well known re- 
sult for pointlike nucleons, (N(k 2} a 2 )\A®\N(kx, ax)) — 

g\uN(k 2 ,a 2 )"/f_ l ^5 :! YU N (kx,ax). We will verify this re- 
lation on quark level, where the axial vector current is 

defined as A®(x) = (u(x) d(x))j t _ l j5 I 2 - (u(x) d(x)) T . To 
operate with matrix elements between either bare proton 
states or bare neutron states we use for the nondiagonal 
cases (a = 1, 2) the assumed isospin symmetry relations, 

cf. H3, {p\u^^ 5 d\n) = (p|u7 Al 7 5 u - d7 M 7 5 d|p) and 

(n|d7^75«|p) = H<^7m75 d - 1*7^75 u\n) [Hi]]. Then, 
taking the solutions of nucleon formula for two-quark op- 
erators, Eqs. (|B5|I . (|B6|I for proton states, and Eqs. (|B7(I . 
(|B8|) for neutron states (see Appendix[BJ) , yields on quark 
level for the bare (isoduplet) nucleon 

(JV(fca >£ r 2 )|l«|JV(Ai j( 7i)> 

= u N (k 2 ,a 2 )j^ 5 —u N (kx,ax) ■ (16) 
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Comparision of l|16[) with the result on effective hadronic 
level yields for the axial charge g v A — \, in fair agreement 
with the value g\ ~ 0.84 deduced from MIT Bag model 
evaluations and neutron /3-decay experiment 62] . 



B. Chiral condensate in nucleon 

The chiral condensate inside the nucleon is related to 
the pion- nucleon sigma term |63| , 

TYl ' - 

a N = 77T7~ (N p hys(k,(j)\uu + dd\N phys (k,a)) , 

ZM TV 

(17) 

where 2m q = m u + mj. A typical value for the pion- 
nucleon sigma term is ctn = 45 MeV [64l IbbT] . The sigma 
term can be decomposed, according to Eq. J2J, into a va- 
lence quark contribution (bare nucleon) and a pion cloud 
contribution (sea quarks and gluons): <tn ~ o% + a N- 
To evaluate a v N we first consider the u quark chiral con- 
densate inside the bare proton. With Eqs. I|B5(1 and 1|B6|) 
one obtains 

{p{k2,(y 2 )\uu\p{k ll (7 1 )) = 2u p (k 2 ,a 2 )u p (k 1 ,a 1 ), 

(p(&2, cr2)\dd\p(ki, (7i)) = 1 u p (k 2 , a 2 ) u p {k\, <n) . 

(18) 

These relations show the momentum and spin depen- 
dence of the chiral condensate inside the bare proton. 
Of course, for a finite-size nucleon there are additional 
momentum dependences for which is accounted for by 
nucleon formfactors. An application of the nucleon for- 
mula to the bare neutron reveals the isospin symmetry 
relations 

(n(k2,(T 2 )\uu\n(ki,ai)) = (p(k 2 , (T 2 )\dd\p(ki, <ti)) , 

{n(k 2 , a 2 )\dd\n(ki, o\)) = {p{k 2 ,a 2 )\uu\p{ki,ai)) . 

(19) 

Accordingly, it is only necessary to compare the findings 
for the proton with results reported in the literature. For 
the special case k\ = k 2 and a\ = a 2 Eq. (|TH|) simplifies 
to 



(p(k,a)\tiu\p{k,a)) = 2 (2.1) , 
(p(k,a)\2d\p(k,a)) = 1 (1.4) . (20) 



2M N 
1 



2M N 

The parenthesized values are the findings of Ref. |66J for 
the valence quark contribution which well agree with our 
results. From (|20|) and isospin symmetry relations i|19|) 
one may now deduce the valence quark contribution to 
the nucleon sigma term within the algebraic approach, 



We compare this result with Ref. [67J , where the valence 
quark contribution to the sigma term has been estimated 
to be <j v n = crjv/(l + Gs fl). By using the given values 
G s = 7.91 GeV~ 2 and f a = 0.393 GeV one obtains a v N = 
20 MeV. Accordingly, our result J2D is, for m q ~ 7 MeV, 
in good numerical agreement with |67| . 

Finally, by assuming that the contribution of the pion 
cloud for the physical proton is the same for the chiral 
u and d quark condensates one can get rid of the term 
ajf by subtracting the chiral d quark from the chiral u 
quark condensate. That means the following approxima- 
tion should be valid 



(p P hys(k,a)\uu - d d\p p hys(k, a)) 
~ (p(k, a)\uu - dd\p(k, a)) = 2M N 



(22) 



a N 



2M N 



^-{N{k,a)\Jlu + dd\N{k,a)) = 3m q . (21) 



where we have used l|18f) . Indeed, the result (1221) is in 

fair agreement with (p p h ys (k, <r)\uu — dd\p p \ vs (k,(7)) = 
2 M N {M s - M E )/m s = 1.3 GeV obtained in 

C. Four-quark condensates 

Four-quark condensates seem to be quite important in 
predicting the properties of light vector mesons within 
the QCD sum rule method 68]. This is related to the 
fact that in leading-order the chiral condensate is nu- 
merically suppressed since it appears in a renormalization 
invariant contribution m q (qq). Therefore, the gluon con- 
densate and four-quark condensates become numerically 
more important. However, the numerical values of four- 
quark condensates are poorly known, and up to now it re- 
mains a challenge to estimate their magnitude in a more 
reliable way. Accordingly, the evaluation of four-quark 
condensates inside the nucleon is an important issue. 
Such quantities have been evaluated in [65| within the 
groundstate saturation approximation, noting the impor- 
tance of four-quark condensates also for properties of the 
nucleon within the QCD sum rule approach. An attempt 
to go beyond the groundstate saturation approximation 
has been presented in |30|, where, by using the Nambu- 
Jona-Lasinio model and including pions and a mesons, 
correction terms have been obtained. Further evalua- 
tions of four-quark condensates beyond the groundstate 
saturation approximation have been performed in Q by 
using a perturbative chiral quark model for describing the 
nucleons. Later, the results of have been used for eval- 
uating nucleon parameter at finite density within QCD 
sum rules [69|. In |9j,|52| lattice evaluations for scalar and 
traceless four-quark operators with non-vanishing twist 
have been reported. In view of these very few results 
obtained so far further insight into such condensates is 
desirable. 

Before considering this important issue we notice a 
general decomposition of four-quark condensates. Let A 
and B two arbitrary two-quark operators. Then the nu- 
cleon expectation value of O — AB can be decomposed 
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as [nJlnl 

(V phys |i£|V phys ) = (i) (V phys |5|V phys ) 
+ (iVph y s|i|A r phys)<S)o + {N phys \AB\N phys ) c , (23) 

where the first two terms refer to the so called factoriza- 
tion approximation, while the last term is a correction 
term to the factorization approximation and describes 
the scattering of a nucleon with B into a nucleon and A, 
i.e. it is a sum over all connected scattering Feynman di- 
agrams TVphys + B — » TVphys + A. The decomposition 
is matched with the decompositions (|24|l . (|25|l and l|45() 
given below, as it is seen in |7l| where we consider an 
explicit example for the vector channel. The first two 
terms in (|23H scale with iV| (N c denotes the number of 
colors, for the moment being taken as a free parameter of 
QCD), while the correction term scales with N c [7^.l73|. 
Dividing both sides of 1)23(1 by N% one recognizes that the 

I 



where q...q is either u...u or d...d (the dots stand for 
r or TA a ). For the chiral condensate we take {qq)o = 
-(0.250 GeV) 3 . The decompositions of Eqs. ^ and 
(|25|l are related to (1 2 ■" > I by means of a Fierz rear- 
rangement; an explicit example for the vector chan- 
nel is given in |71|. The last term on the rig ht side 
of Eqs. (|24|l and (|25|l is a correction term |30( to the 
groundstate saturation approximation |65| . describing 
the scattering process 7V p hy S + q...q — ► N p \ lys + q...q. 
To get an idea about the magnitude of these correction 
terms we consider two typical examples. The factor- 
ization approximation l|24[) (i.e. without the correction 
term) yields for the scalar channel ( V p hy S | qq\N p hy S ) = 
—0.173 GeV 4 . In [3(| the correction term to this ground- 
state saturation approximation has been found to be 
(-^phys | qq qq I Aphys) C = 0.011 GeV 4 . As another example 
we consider the vector channel with Gell-Mann matrices. 
The factorization approximation l|25|l (i.e. without the 
correction term) yields (Np^q-y^q ^^\ a q\ V p hys ) = 
0.335 GeV 4 , while the correction term in [3(| is 



last term has to be considered as a correction term of the 
order l/N c jzJEl- 

That means that a factorization of 
four-quark operators in a cold medium is consistent with 
the large- iV limit [72LI73T ]. a statement which is also valid 
in vacuum l74l. 



1. Flavor-unmixed four-quark condensates 

We start our investigation with the flavor unmixed 
four-quark condensates and consider the general expres- 
sion of two different kinds of flavor unmixed condensates 
inside the nucleon, namely condensates without and with 
Gell-Mann matrices A a (a = 1, 8 are the Gell-Mann in- 
dices, which should not be confused with the color indices 
(in roman style) a, b,c (later also i...,n) of quark fields) 

[Hlfllz! ' 



i 

(V phys |f7 Al A a gf 7 ^A Q ( 7|V phys ) c = -0.139 GeV 4 . Accord- 
ingly, the correction to the groundstate saturation ap- 
proximation in the scalar channel turns out to be less 
than 10 percent, while in the vector channel with Gell- 
Mann matrices it is about 30 percent. As we will see, 
from Q24H and l|25|l the valence quark contribution can 
be extracted in a unique way. 



Now we evaluate the valence quark contribution of 
four-quark condensates, and start to consider the u quark 
inside the bare proton. Application of the nucleon for- 
mula ||HJ) with the composite proton field operator l|12|) 



(Nphyslglig gTagliVphys) = - 



Tr(r 1 )Tr(r 2 )--Tr(r 1 r 2 ) 



(qq)o {N pYlys \qq\N phys } 



1 

16 



Tr(ri)Tr( 7 T 2 ) + Tr(r 2 )Tr( 7 ' i r 1 ) - ~Tr(r l7 T 2 ) - ^Tr^^) 
X (q~q)o (N p hys\qJnq\N phys ) + (N phys \qT 1 qqT 2 q\N phys ) c , 



(24) 



„ |f riA a gtr 2 A a (7| V p hys) = --Tr(r!r 2 ) 0q) o (N phys fiq\N phys ) 



Tr(r l7 ^r 2 ) + Tr(r 2 rr 1 ; 



(fg) (V phys |f 7 ^g|V phys )+ (V phys |fr 1 A a gfr 2 A a g|Af phys ) c , (25) 
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yields 



{p{k 2 ,a 2 )\uT 1 u uT 2 u\p(k-L,ai)) 

= U^ 2 (fc 2 ,cr 2 ) (7o)/3 2 a 2 (7o)ai/3i ^(fajITi) 

x (Ti) a0 {T 2 y 5 J d 3 n e ikiri / d 3 r 2 e"* 21 " 2 



V£ 2 (r 2 ,0) 



« a "Js S ' ( ri '°) 



By inserting the expression given in Eq. (|B9|) in the Ap- 
pendix El into H2fijl one obtains for the bare proton 



|0) . 

(26) 



(p{k 2 ,(72)\uT 1 u uT 2 u\p(k 1 ,ai)) = ^ (uu)o ^3 Tr(ri)u p (fc 2 , (J 2 )T 2 u p {ki, cr x ) + 3 Tr(r 2 )u p (/c 2 , (J 2 )Tiu p {k u a{) 



-M p (fc 2 ,cr 2 )r 1 r 2 Mp(fci,(T 1 ) - w p (fc 2 ,cr 2 )r 2 r 1 Up(fc 1 ,cr 1 ) 



(27) 



r 



In a similar way one obtains for the four-quark conden- 
sates involving Gell-Mann matrices 

(p(*a, CT 2 )|M\A a w tr 2 A a {t|p(fci, cri)) 
= -^-(Jlu)o\u p (k 2 ,a 2 )T 1 r 2 Up(k 1 ,(T 1 ) 



+u p {k 2 ,a 2 )Y 2 T 1 u p (k u cr 1 ) 



(28) 



For the d flavor we have 



{p(k 2 ,a 2 )\dTid dT 2 d\p{kx,(T\)) 

= ^ (P( fc 2,0" 2 )|uriu ur 2 u|p(fci,cri)) , 



(29) 



{ P (k2,a 2 )\dT 1 X a d dr 2 A a d|p(fci, ai)) 

= g (p(fc 2 , ( T 2 )|^r 1 A Q u^r 2 A Q u|p(fci,ffi)) • (30) 

An analog evaluation of these four-quark operators inside 
the bare neutron gives 

{n(k 2 , a 2 )\uTiu uT 2 u\n{ki,oi)) 
= {p{k 2 ,a 2 ){dT l dlT 2 d\p{k ll G 1 )) , (31) 

{n{k 2 , er 2 )|uriA a u uT 2 X a u\n(ki, cri)) 
= (p(k 2 ,a 2 )0T 1 X a dlr 2 X a d\p(k 1 ,a 1 )) , (32) 

which, like in the case of chiral condensate, reflect the 
isospin symmetry. By interchanging u <-» d on both sides 
in Eq. I|31|l and 11321) one also gets the d flavor four-quark 
condensates inside neutron. The results Ij27(l - Ij32(l for 
the four-quark condensates inside the bare nucleon distin- 
guish between proton and neutron, and they also contain 



the dependence of these condensates on the momentum 
of the (pointlike) nucleon. Of course, as in the case of 
two-quark matrix elements, for a finite-size nucleon there 
are additional momentum dependences implemented in 
formf actors. 

We compare now these findings of Eqs. (|27|l - i|32|) with 
Ref. i- e - we se t fa — k 2 and o\ = a 2 , and average 
over proton and neutron to get the nucleon condensates, 
2(N\0\N) = (p\6\p) + (n\0\n). First we consider the 
case of scalar four-quark condensate, i.e. Ti = T 2 = 1. 
Then, for the bare nucleon one obtains from Eqs. 127|l . 
(PI and (El) 

(N(k,a)^q\N(k,a)) = y (f$> M N , (33) 

while, according to the first term in Eq. I|24|l (the sec- 
ond term vanishes in this special case), for the physical 
nucleon the result 



(N p h ys (k,(j)\qq qq\N phys (k,a)) 



11 (a v N air 

— (qq)o M N I — + — 
u \ m q m q 



— - (qq) M N — 



(34) 



has been obtained in [6^|. For the last line of Eq. I|34l) 
we have used the decomposition <7n — crjj + ajf, which 
has already been considered in section IIVBI Compar- 
ing both results, Eqs. I|33[l and l)34|l. one recognizes, by 
means of relation <7 v N /m q = 3 (cf. (63 and Eq. J2U), 
that the result (|3^1) is nothing else but just the valence 
quark contribution of the scalar four-quark condensate 
inside the nucleon; it is in agreement with the separated 
valence quark term of Eq. (|34f) • 

Due to its importance and its instructive property we 
will also consider the case Y\ = l,r 2 = j p . From 
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Eqs. (EH, (EH and gTJ we find 

(7V(fc,a)|fgf 7p (7|iV(fc, ( T)) 

5 

= j{M)o (j)^ p u N (k, a) , (35) 

which is the contribution of the bare nucleon, i.e. the va- 
lence quark contribution. To compare it with the corre- 
sponding result of Ref. J6f| we first deduce from Eq. I|24() 
that 

(N phys (k, a)\qq qj p q\N phys (k, a)) 

5 

= g(OT)o (Np h y S (k,a)\q-f p q\N ph y S (k,(7)) . (36) 

From H36J) we have to extract the valence quark contri- 
bution by virtue of Eq. @ (with ~ 1) 

(N phys (k, cr )l<77p'7|AW s (fc, er)) 

= (N{k,(r)\h P q\N(k,(T)}+J2€ (Nmr\iijpq\Nmr) . 

n 

(37) 

The first term on the right side, which is in fact an aver- 
age over proton and neutron, is the valence quark term 
we are interested in, while the second term is the pion 
cloud contribution. With isospin invariance one obtains 

(N{k,a)\h P q\N{k,a)) 

= ^(P(k,<j)\uj p u\p{k,a)) + {p{k,a)\d'ypd\p{k,a)) S j 

= \{^ "p( fc ' a )lp u p{ k i cr ) + 1 U p (k, <j)j p u p (k, cr)^ 
3 

= -u N {k,a)^ p u N {k,a) , (38) 

where for the last term we have set u p — un because of 
M p = M N . By inserting into and using we 
obtain 

(N ph y S (k, a)\qq qj p q\N phys (k, a)) 
5 

= -{qq)aUN(k,a)jpU N (k,a) 

+ lm)oY, ( l ) n( Nn7T \H<i\ Nn7T ) ■ ( 3Q ) 

n 

The first term on the right side of Eq. (|39|l agrees with 
our result l|35|l . while the second term on the right side 
of Eq. H39|) is a factorization approximation of the ex- 
pression Y^, <t>n {N ni:\qqq^ p q\N nn) . From that it becomes 

n 

obvious that our result (|35|l is nothing else but just the 
valence quark contribution of 13911 . 

Other combinations of Clifford matrices, like Ti — 75 
and T2 = 7 P , with or without Gell-Mann matrices, can 



be evaluated and compared in the same way |76j. This 
means that within the algebraic approach (JHJ) for evaluat- 
ing bare nucleon matrix elements we find an agreement 
for all of the flavor-unmixed four-quark condensates if 
one takes from the corresponding results of Ref. [65| the 
valence quark contribution, for instance by means of the 
decomposition on = a N + a N . Therefore, one actually 
may consider our evaluation as a re-evaluation of the va- 
lence quark contribution of the four-quark condensates of 
Ref. |£5j and a confirmation of their results within an in- 
dependent microscopic approach (quark-diquark picture 
of the nucleon). But we have to be aware that such 
an agreement between our algebraic approach and the 
factorization approximation applies only for the valence 
quark contribution of nucleon matrix elements. Espe- 
cially, such an agreement is not expected when taking 
into account the pion cloud contributions according to 
Eq. 

Having found the remarkable agreement with valence 
quark contribution of the factorization approximation 
it becomes interesting to compare our results also with 
other evaluations in the literature. However, it turns 
out that a comparison with the recent lattice data of 
Ref. 91 is quite involved since in [[j the condensates 
have been evaluated at a renormalization scale of about 
lattice — ^ GeV 2 , which is considerably higher than our 
renormalization point of about /j? ~ 1 GeV 2 (our renor- 
malization point is hidden in the chiral condensate, i.e. 
(qq)o(^ 2 ))- To scale the lattice data from 5 GeV 2 down 
to the hadronic scale of 1 GeV 2 one needs to know the 
matrix of anomalous dimension for all of the four-quark 
operators which accounts for the effect of operator mixing 
among the four-quark condensates. This operator mixing 
may change considerably the numerical values and even 
the signs of the four-quark condensates given in [l| . An- 
other method which seems also capable to compare our 
results with lattice data would be to scale our renormal- 
ization point n 2 up to the lattice scale At 2 atticc - Such a 
procedure requires, however, the knowledge of renormal- 
ization scale dependence of the chiral condensate. Allto- 
gether, performing these procedures is out of the scope 
of the present paper and we therefore abandon a com- 
parison of our results with the ones given in Ref. 0. 

In view of the mentioned points and especially in view 
of an acceptable lucidity of our paper, we restrict our- 
selves to a comparison with the recently obtained four- 
quark condensates of Ref. ||| . Due to the specific nota- 
tion for the four-quark condensates choosen in Ref. Q we 
list our results for the valence quark contribution of scalar 
four-quark condensates for all channels in the same way 
as in Ref. 8] . Our results for the valence quark contribu- 
tion can be obtained from Eqs. lt77|) - by averaging 
over proton and neutron (k\ = &2, o\ = 02): 
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{N\\im-\^ a qtX a q\N) 

2. ± 1~ - 
{N\ - qj 5 q 9759- 2975A a gg7 5 A a g|7V) 

(N\~hnqh»q ~ ^7^ a ^7 M A a (z|iV) 

2 1 

(^1 3 97^759 qi^lM ~ ^1 »l$X a qq^ q\N) 

(Nl-q-a^q^q - ^a^^q^ X a q\N) 



-0.0733 GeV 4 (-0.117 GeV 4 ) , 
-0.0147 GeV 4 (-0.0567 GeV 4 ) , 
-0.0586 GeV 4 (-0.0582 GeV 4 ) , 
+0.0586 GeV 4 (+0.0567 GeV 4 ) , 
-0.176 GeV 4 (-0.356 GeV 4 ) . 



(40) 
(41) 
(42) 
(43) 
(44) 



r 



The values parenthesized are the results for these conden- 
sates as given in Ref. |8(, but for the physical nucleon, 
i.e. for a nucleon which contains the valence quark, sea 
quark and gluon contributions [Tt]]- Since we have com- 
pared our evaluated valence quark contribution with the 
total contribution for the physical nucleon of Ref. |8| it 
becomes obvious that one may actually not expect a per- 
fect numerical agreement. The more interesting fact is 
to notice that the valence quark contribution for the vec- 
tor and axial vector channel, (|42|l and 143fl . respectively, 
agrees very well with the total contribution for the phys- 
ical nucleon. For the scalar, axial scalar and tensor chan- 
nel the signs for the valence quark contribution and total 
contribution of Ref. |]| agree, while the numerical values 
differ. This illustrates that the sea quark and gluon con- 
tributions are expected to give noticeable contributions. 



2. Flavor-mixed four-quark condensates 

For the flavor-mixed four-quark condensates the gen- 
eral decomposition reads [65| 



the vector and axial vector channel with Gell-Mann ma- 
trices). 

Now we consider the valence quark contribution of the 
flavor mixed four-quark condensates, i.e. the contribu- 
tion of the bare nucleon. Application of our nucleon for- 
mula JSJl with (|12fl yields for the bare proton 

{p{k 2 ,a 2 )\uT 1 u dT 2 d\p(k ll ai)) 

= - 1 2 (uu) Tr(ri) u p (k 2 , a 2 ) T 2 u p (ki,ai) 



1 (dd) Tr(r 2 ) u p (k 2 ,a 2 )Yiu p {ki,<Ti) 



(46) 



{p{k 2 , a 2 )\llTi\ a u dT 2 \ a d\p{k u ai)) = , (47) 
while for the bare neutron we find 

(n(fc 2 , cr 2 )|?iriu dr 2 d|n(/ci, (Ti)) 

= - ( 1 (dd) Tr(ri) u n (k 2 , a 2 )T 2 u n (k 1 ,a 1 ) 



(N phys \uTiu dT 2 d\N phys ) 

= o(qq)o (N phys \qq\N phys ) Tr(r x ) Tr (r 2 



+ (qq)o {Nphyslqj^Nphys) 



Tr(ri) Tr( 7 T 2 ) 



-Tr(r 2 ) Tr( 7 ^r 1 ) 



{NphysluTxu dT 2 d\N phys ) c , 

(45) 



where the isospin symmetry relations l|19Jl have been 
used. The last term in (|45|l is a correction term to the fac- 
torization approximation, describing the scattering pro- 
cess JVphys + dT 2 d — -> iVphys + uTiu. The decomposition 
(|4*5jl is, like (H3J) and J25J|, matched with Eq. |23I> by 
means of a Fierz transformation. The flavor-mixed con- 
densates with Gell-Mann matrices vanish in the factoriza- 
tion approximation, (N p Y 1 y S \uT 1 \ a u dT 2 X a d\ N p h ys ) = 
[65| (in [7^| we have found small corrections to the fac- 
torization approximation of flavor-mixed condensates for 



■ 2 (uu)o Tr(r 2 ) u n {k 2l <7 2 ) T x u n {ki,a{] 



(48) 



(n(fc 2 ,(T 2 )|^riA a wdr 2 A a (l|n(fc 1 ,cr 1 )> = 0. (49) 

The Eqs. I)46|) and (|48|l are generalized expressions of the 
factorization approximation given in |65j | because of the 
dependence on nucleon momentum and the distinction 
between proton and neutron. The results of Eqs. (|47|l and 
(|49|l are in agreement with the factorization approxima- 
tion |65|| . but, as mentioned above, beyond the factoriza- 
tion approximation these condensates are nonvanishing 

M 

To compare the obtained results of Eqs. PH)) and 
with the factorization approximation, i.e. neglecting the 
correction term in l|45|) , we consider the special case k\ — 
k 2 , <7i = a 2 and Ti = T 2 = 1. For the bare nucleon 
one obtains by averaging over the bare proton and bare 
neutron, 



{N(k,cr)\Tiu dd\N{k,cr)) = 6 (fg) M N 



(50) 
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This result has to be compared with the corresponding 
finding of j6^ which according to Eq. 145f) reads 



(N phys (k,a)\uu dd\N phys (k,a)) = 2 (qq) M N 



UN 



2 (qq)o M N ^ + 

Tfln Win 



(51) 



where in the last expression we have used the decomposi- 
tion crjv = a^+ajf. By means of the relation a v N jm q — 3 
(cf. [Ell and E q- <ED) the result ^ is in agreement with 
the separated valence quark contribution of l|51|) . As in 
the cases considered in the previous subsection IIV C II 
such an agreement with Ref. (6j| can be achieved for 
all the other combinations of Fi and Y 2 of the Clifford 
algebra. 

As in the case of flavor-unmixed four-quark operators 
we compare our findings for the valence quark contribu- 
tion of flavor- mixed condensates, (|46[) and (|48H . with the 
total result for the physical nucleon of Ref. y| to exam- 
ine the magnitude and sign of our results. According to 
Eqs. (|46|l - l|49|) only the valence quark contribution of the 
scalar channel does not vanish. Its numerical magnitude 



{N\\mdd- Ul\ a ud\ a d\N) 
-0.0586 GeV 4 (-0.094 GeV 4 ) , 



(52) 



turns out to be comparable with the evaluation of [8| for 
the total contribution of the physical nucleon given in the 
parentheses in (|52|1 . The numerical difference in magni- 
tude is caused by sea quark and gluon contributions. 

To summarize this section, we have evaluated the va- 
lence quark contribution to flavor-unmixed and flavor- 
mixed four-quark condensates for the u and d flavor 
inside proton and neutron. The results for the flavor- 
unmixed operators are given by the Eqs. I|27|) - (|32() . and 
the results for the flavor-mixed operators are given by 
the Eqs. (|4*uT) - (|4*9"|) . We have seen that our findings 
for the four-quark condensates within the algebraic ap- 
proach, which is by far a different method than the used 
one of Ref. [6{| , are in agreement with the large- N c limit 
[l2>lz3 and with the results of Ref. J6{| when taking from 
there the valence quark contribution only. It seems ad- 
missible, especially in view of the agreement with valence 
quark contribution of factorization, that the nucleon for- 
mula yields reliable results for the valence quark contri- 
bution of four-quark condensates inside the nucleon. 



V. SIX-QUARK CONDENSATES 

Six-quark condensates become important mainly for 
two reasons. First, in the operator product expansion 



(OPE) of current correlators one usually takes into ac- 
count all terms up to the order of the four-quark conden- 
sates and neglects the contributions of higher order. Such 
an approximation may work or may not work, depending 
on the specific physical system under consideration. Ac- 
cordingly, one has to be aware about the contribution of 
the next order to decide how good such an approximation 
is. This is also necessary for the more involved case of 
finite density, where a Gibbs average over all hadronic 
states of the correlator under consideration has to be 
taken. Indeed, a very recent estimate of such higher con- 
tributions beyond the four-quark condensate for the nu- 
cleon correlator in matter underlines also the importance 
of an estimate for the six-quark condensates inside the 
nucleon |79| . And second, it is well known that instan- 
tons give rise to corrections to the Wilson coefficients of 
six-quark condensates [80Tl8l| . These corrections cause a 
substantial enhancement of the vacuum contribution of 
six-quark condensates in the OPE of current correlators. 
To investigate such current correlators at finite density 
implies the knowledge of the nucleon matrix elements of 
six-quark condensates. Here, after getting confidence on 
our proposed approach in the previous sections, we will 
use the nucleon formula to evaluate the valence quark 
contribution of six-quark condensates inside the nucleon. 

Within our algebraic approach by using the nucleon 
formula © with i|12|) we obtain for the u flavor inside 
the bare proton 

{p(k 2 ,a 2 )\uTi u uT 2 u uT 3 u\p(ki,ai)) 
= u^ 2 (k 2 ,a 2 ) (70 W 2 (7o)ai/3! u^{ki,ax) 



x (T,)^ (T 2 y S (r 3 ) eC | d 3 ne lkiri J d 3 r 2 



-ik 2 r 2 



<0| 



V£ 2 ( r 2>o), 



u a u-p u y ?4 «r > % ( r i>°) 



|0) . 

(53) 



The commutator-anticommutator is given in Eq. (|B11|) 
in Appendix iBlfor a more general case. According to this 
result the six-quark condensate inside the bare proton is 
reduced to a four-quark condensate in vacuum. We note 
one of these four-quark condensates in vacuum saturation 
approximation Q 

{u% w 7 Ug U e )o = 

x (sfrSs, V k 8 lm - 5f3 e 5 lS 8 im S kl ) . (54) 



When evaluating all of the four-quark condensates of 
Eq. i|Bllfl in the same way one obtains, by using the 
normalization (|13fl . 
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(j)(k 2 ,a 2 )\uTiuuT 2 uuT 3 u\p(ki,ai)) = (' 



Up(fc 2 , cr 2 ) (ri r 2 r 3 + r 3 r 2 + r 2 r x r 3 ) up(fci, 01) + u p (k 2 ,a 2 ) (r 2 r 3 ri + r 3 r x r 2 + r 3 r 2 r x ) 01) 

-3u p (k 2 ,a 2 )T 1 T 3 u p (k 1 ,ai) Tr(r 2 ) - 3u p (k 2 ,(X 2 )T 2 T 3 u p (ki,(Xi) Tr(r x ) - , Zu p {k 2 ,cr 2 )Y 1 T 2 u p (k- i _,cri) Tr(r 3 ) 
-3u p (fc 2 ,cr 2 )r 3 r 2 Up(fci,cri) Tr (ri) - 3u p (k 2 ,a 2 ) (r 2 Ti) u p (ki,<7i) Tr (T 3 ) 
-3Mp(fc 2 ,cr 2 ) (r 3 ri)up(fci,cri) Tr(r 2 ) + 9lIp(fc 2 ,c7 2 )r 3 Up(fci,£ri) Tr(ri) Tr(r 2 ) 
+9up(fc 2 ,cr 2 )r 2 Mp(A;i,o'i) Tr(r x ) Tr(T 3 ) + 9 u p {k 2 , a 2 )T 1 u p {k 1 , a\) Tr (r 2 ) Tr (r 3 ) 

-3u p (k 2 ,a 2 )T 3 Up(ki,ai) Tr (r x r 2 ) - 3 u p (k 2 , a 2 )T 2 u p (hi, cr x ) Tr^xTa) - 3 Up(fc 2 , (7 2 )riu p (fci, a x ) Tr (r 2 r 3 ) 



(55) 



r 



For the d flavor six-quark condensate we get 

(p(k 2 , cr 2 )|rfri d dT 2 d 3r 3 d\p(ki,ai)) 
= ^(p(k 2 ,(T 2 )\uTiuur 2 uur 3 u\p(ki,cri)) . (56) 

These findings are, to the best of our knowledge, the first 
attempts to estimate the size of six-quark condensates 
inside a (bare) nucleoli. An analog evaluation for the 
bare neutron reveals the isospin symmetry relations 

(n(k 2 ,a 2 )\uTi uuT 2 uuT 3 u\n(ki, oi)) 

= (p(k 2 ,a 2 )\dTiddT 2 ddT 3 d\p(ki,(7i)) , (57) 

(n(k 2 , a 2 )\dT\ d dT 2 d dT 3 d\n{k\, ci)) 
= {p{k 2 ,(x 2 )\uT\u uT 2 uuT 3 u\p(ki,cri)) . (58) 



Finally, by averaging over the proton matrix elements, 
Eqs. (|55(l and (|56|) . and the neutron matrix elements, 
Eqs. ij57)l and one gets the six-quark condensates 

inside a bare nucleon. For instance, the six-quark con- 
densate for the scalar channel is found to be 



55 

(N(k,a)\t qtqqqWk,*)) = y(f<?)o M N . (59) 



Further, we present results for six-quark condensates 
which contain Gell-Mann matrices. Note that only nu- 
cleon matrix elements of colorless operators do not van- 
ish, i.e. only two Gell-Mann matrices may occur. With 
the general result i|Bllfl in Appendix [B] and normaliza- 
tion (|13|) we obtain 



(p{k 2 ,a 2 )\uTi\ a u uT 2 \ a u uT 3 u\p{k 1 ,a 1 )) = j—^ (uu)l 

^Up(fc 2 ,er 2 ) (riT 2 r 3 + Tir 3 r 2 + r 2 rir 3 )wp(fci,o-i) + -^-u p (k 2 ,a 2 ) (r 2 r 3 ri + r 3 rir 2 + r 3 r 2 r :L )'Up(/ci,o-i) 

-16Mp(fc 2 ,cr 2 )r 2 riUp(fci,cri) Tr(r 3 ) - 16up(fc 2 ,cr 2 )rir 2 ?Zp(£;i,c7i) Tr(r 3 ) - 16wp(fc 2 ,(7 2 )r 3 Up(fci,cri) Tr(rir 2 ) 



(60) 



r 



For the d flavor we get 



(p(k 2 , CT 2 )|dTiA a d dT 2 X a ddT 3 d\p{k u oi)) 
(p(k 2 ,a 2 )\fir 1 X a u%T 2 X a uTiT 3 u\p(ki,(Ji)) . (61) 



Finally, evaluating these operators inside the bare neu- 



tron we obtain the isospin symmetry relations 

(n(k 2 , a 2 )\uTi\ a u ?ir 2 A a u uT 3 u\n(ki, cri)) 
= ( P {k 2 ,a 2 )\dT 1 X a dlT 2 X a d3r 3 d|p(fci,(7i)) , (62) 
(n(k 2 , c7 2 )\2r 1 X a ddT 2 X a ddY 3 d\n{k x ,(Ji)) 

= (p(fc2,CT 2 )|^riA a u^r 2 A Q u^r 3 u|p(fci,(Ti)} . (63) 

As before, by averaging over the proton matrix elements, 
Eqs. (|60|l and (|61|) . and the neutron matrix elements, 
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Eqs. H62|) and llG.'ill . one obtains the six-quark conden- 
sates containing Gell-Mann matrices inside a bare nu- 
cleon. The presented findings for six-quark condensates 
inside the bare nucleon provide basic results for further 
investigations beyond the order of four-quark conden- 
sates within the QCD sum rule approach for the nucleon 
in matter (for the nucleon sum rule in vacuum see 
and for the nucleon sum rule in matter see |64l l65l l82fl ) . 
Investigations aiming at predictions beyond the order 
of four-quark condensates, however, imply in addition 
to the evaluation of the six-quark condensates also the 
knowledge of the Wilson coefficients for all of these six- 
quark condensates. So far, these coefficients in the OPE 
for the nucleon correlator have been determined up to 
the order of the four-quark condensates [zi| • 



VI. SUMMARY 

An algebraic approach for evaluating bare nucleon ma- 
trix elements has been presented. The supposed nucleon 
formula © relates bare nucleon matrix elements to vac- 
uum matrix elements and, therefore, new parameters for 
evaluating them are not needed. A feature of the alge- 
braic method is that the valence quark contribution of 
two-quark, four-quark and six-quark condensates can be 
evaluated on the same footing. One aim of the present 
paper is to demonstrate how the nucleon formula works 
and to test it in several cases. In doing so, the nucleon 
has been considered as a composite pointlike object, de- 
scribed by a valence quark and a valence diquark ap- 
proximated by two classical Dirac spinors. Neither sea 
quarks nor gluons, or in a hadronic language no meson 
cloud, have been taken into account here. Accordingly, 
the results presented have to be considered as pure va- 
lence quark contribution to the matrix elements under 
consideration. 

A consideration of the electromagnetic current i|15|) 
and the axial vector current (|16fl for the bare nucleon 
reveals the expected current structure for a pointlike ob- 
ject. We have evaluated then the valence quark contri- 
bution of the chiral condensate (|21[l , finding the relation 
(t v n = 3 m q which turns out to be in numerical agreement 
with results obtained in an earlier work [67|. 

Furthermore, the valence quark contribution of four- 
quark condensates has been investigated. Our results 
are given in Eqs. (|27() - (|32|l for flavor- unmixed opera- 
tors, and in Eqs. (|4(j|l - l|49|> for flavor-mixed operators. 
In the special case fci = k2, <J\ = 02 we find an interest- 
ing agreement with the groundstate saturation approxi- 
mation explored in Ref. j65| if one separates the valence 
quark contribution of four-quark condensates from that 
results. In this respect our approach yields an indepen- 
dent re-evaluation and confirmation of the results of Ref. 
[65| . because both methods are different from the con- 
ceptional point of view, which is already interesting in 
itself. Even more, our algebraic approach presented re- 
covers the dependence of condensates on momentum for 



a pointlike nucleon and distinguishes between proton and 
neutron matrix elements. In this respect it goes beyond 
the common factorization approximation. In this con- 
text it is worth to underline that the agreement between 
our algebraic approach and the groundstate saturation 
approximation has been found for the bare nucleon, and 
not for the physical nucleon. In Eqs. (|40|) - l|44|) and 
(|52|l we have compared our results with values of four- 
quark condensates inside the physical nucleons recently 
obtained within a chiral quark model 

As a further application of nucleon formula we have 
presented results for six-quark condensates inside the 
bare nucleon, given in Eqs. - and Eqs. (|oT])l - 
(|63|l . respectively. These values obtained are, to the best 
of our knowledge, the first evaluation of six-quark con- 
densates inside (bare) nucleons. Finally, in Eqs. i|21[l . 
(|33|l and (|59|l we have given more explicit examples for 
the scalar channel of two-quark, four-quark and six-quark 
operators, respectively, inside the bare nucleon, showing 
up an interesting alternative change in the algebraic sign 
from two-quark to four-quark and from four-quark to six- 
quark condensates. 

A remark should also be in order about nucleon matrix 
elements of gluonic operators. Evaluating such operators 
within the algebraic approach presented requires, in gen- 
eral, the implementation of gluonic degrees of freedom 
into the composite nucleon field operator H12|l . Such an 
implementation might be provided by the quark-gluon 
interpolating nucleon field discussed in another context 
in H|. 

The algebraic approach can be extended into several 
directions. First, the description of the proton core with 
the field operator (|12f> . and the corresponding one for the 
bare neutron, can be improved, e.g. by implementing an 
effective potential between the valence quarks. And sec- 
ond, the pion cloud of nucleon, accounting for virtual 
sea quarks and gluons inside the physical nucleon, can 
be implemented within the Tamm-Dancoff method. To 
get an algebraic approach also for such a case one has to 
combine the nucleon formula with the soft pion theo- 
rem (IA8II . This implies the evaluation of the coefficients 
4>n in within the renormalizable pion-nucleon inter- 
action Hamiltonian, which is therefore a topic of further 
investigations. Accordingly, for the time being the ap- 
plication of nucleon formula (JSJ in combination with the 
field operator (|12fl has to be considered as a first step in 
determining more accurately nucleon matrix elements of 
quark operators. In summary, we arrive at the conclusion 
that a reliable evaluation of quark operators inside nu- 
cleons can be based on a purely algebraic approach. This 
triggers the hope that predictions of in-medium proper- 
ties of hadrons become more precise in future. 
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APPENDIX A: SOFT PION THEOREM 

In this Appendix we recall a soft pion theorem relevant 
for our purposes to show the similarity in derivation and 
final form of it with the nucleon formula (J8J. Let us 
consider the general pion matrix element of an operator 
(x) which in general may depend on space and time 



(7r b (p 2 )\d(x)\7T a ( Pl )) 

= ^Z- 1/2 J d 4 x ie -^ (D xi +ml) 

x (n b (p 2 )\ T W (6{x) a ( Xl )) |0> , (Al) 

where the LSZ reduction formalism has been applied on 
pion state \-K a (pi)). Here, the letters a,b — 1,2,3 de- 
note isospin indices. The normalization of pion state 
is (n b (p 2 ) \ir a ( Pl )) = 2£ pi (2 7 r) 3 (5( 3 )(p 1 - p 2 ) S ab , where 
E pi = y/pl + rri*. The wave function renormalization 

— 1/2 

constant is < Z v < 1. The normalization of nonper- 
turbative QCD vacuum is (0|0) = 1. Here, x\ — (ri,ii) 
is the space-time four-vector, and Ty/ denotes the Wick 
time ordering. The states |7r a (pi)) are, of course, on-shell 
states, i.e. solutions of the Klein-Gordon equation for 
nonintcracting pions, while the field operator (p a , in gen- 
eral, is the interacting field, i.e. it is off-shell. The four 
momenta in (|A1(I are on-shell, i.e. p\ = p\ = m 2 . The 
soft pion theorem is valid for a noninteracting pion field 
(i.e. Z v ' = 1), wich is a solution of the Klein-Gordon 
equation [D x + m^) (p a {x). In order to be complete in 
the representation we will also specify the noninteracting 
pion field operator 



(A2) 



where the creation and annihilation operators obey the 
following commutator relations 



a a ( Pl ), a b \ P2 ) 



b a {vi),V>\p 2 



2E P1 (2w) 3 S^( Pl -p 2 )S ab . (A3) 



Accordingly, |tt q (p)) = a a ^(p)\0). From |(X2} and $M 
we deduce the equal-time commutator for the nonintcr 
acting pion fields, 



0°(n,t) , do<p b (r 2 ,t) 



= i<5^(n -r 2 ) S ab . (A4) 



In addition, the soft pion theorem is only valid in case of 
vanishing four- momentum p^ — > which implies = 0. 
Then we get 

km {n b (p 2 )\d(x)\7r a ( Pl )) 

Pl^O 

d 4 Xl D X1 (ir b (p 2 )\ T w (6(x) a (xi)) |0) 



|0) S(t-t!) 



(A5) 



In the last line we have used the equation of motion for 
the massless pion field, i.e. Xl (p a (xi) = 0. Now, the 
PCAC hypothesis, which asserts a relation between the 
axial current and the pion field (/„. ~ 92.4 MeV is the 
pion decay constant), 



A°Jx) = -Ud^ a (x) , 



(A6) 



is inserted into Eq. (|A5|) (by means of the field equation 
for the noninteracting pion it becomes obvious from (|A6() 
that in the limit of vanishing pion mass the PCAC goes 
over to a conserved axial vector current (see also 59])). 
The axial vector current A"(x) obeys the well known 
current algebra commutation relations 12] which directly 
leads to the soft pion theorem relevant for our purposes 

M 



lim {n b {p 2 )\0(x)\ir a ( Pl )) = - / d A Xl 



X (7T b (p 2 )\ 



6{x) , A a ( Xl ) 



|0) Sit -h 



(A7) 



One may apply the same steps as before on the other 
pion state as well, ending up with the soft pion theorem 

lim lim (7r 6 (p 2 )|6(a;)|7r a (p 1 )) 

P2—"0 Pi—*0 

= J2 j ^1 j d4x 2 S(t - h) 5{t - t 2 ) 



(0| 

1 

72' 

J 7T 



A b (x 2 ) 



O(x) , A a { Xl ) 



10) 



0{x) , Q\ 



|0). 



(A8) 



where Q\ is the (time independent) axial charge, Q a A — 
J d 3 r Ag(r, t). It seems expedient to emphasize that the 
QCD quark degrees of freedom were not necessary for 
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deriving the soft pion theorem. If one expresses the axial 

vector current by quark fields, A^(x) — ^(r a /2)7 m75 ^ 

with ^> = (u d) T , and the pion fields as well by means 
of their interpolating fields (i.e. composite quark fields 
which have the quantum numbers of pions), then the re- 
lation l|A6() . the current algebra and, therefore, the soft 
pion theorem (|A8J| can also be established within QCD 
degrees of freedom. This theorem can then also be used 
for evaluating pion matrix elements of quark operators 
(see, for instance, 0,122]). Summarizing, the soft pion 
theorem is valid for a noninteracting pion field with van- 
ishing pion four-momentum. In many applications such a 
restriction is not problematic since the pion mass is small 
compared to a typical hadronic scale of about 1 GeV. 



APPENDIX B: EVALUATING EQ. lt53l) 

To evaluate Eq. I|53[) we start with the case of two 
quark field operators. The nucleon formula (JSJ with H12|) 
yields 

(p(k 2 ,a 2 )\u a U^|p(fci,(Ti)) 

= u^(k 2 ,a 2 ) (7o)/3 2Q2 (70)01/81 u^(ki,ai) 

x f d 3 n f d 3 r 2 e ikiri e- 4k2r2 



|0) .(Bl) 



Inserting the proton field operator 112(1 and using 
[AB,C]- = A[B,C] + - [A,C]+B for the commutator 
we obtain 



u a u ,tp p (ri,0) 



j^* a b c ( —a T 



u al (ri)C 7 5d (ri) 



x S ic (70 W <5 (3) (ri)C ■ (B2) 



In the same way we get 



-I + 



x 6 CC ' J (3) (ri) <5 (3) (r 2 ) 

a ' T (r 2 ) C75 d h '(r 2 )) (V T (n) C75 rf b (rx)) 



(B3) 



Integrating over both delta-functions and then using the 
normalization (|13[1 we obtain 



d 3 n e lkiri / d 3 r 2 e^ lk2r2 



x (0| 



^p 2 (r 2 ,0), 



u a Up,4> p (ri,0) 



|0) 



= 2 (7o)/3ai (jo)aa 2 , 

and with HBH 



(B4) 



(p{k 2l a 2 )\u a Up\p(ki,ai)) = 2 Up(k 2l a 2 ) u^{k\,a\) 



(B5) 



We note an analog relation for the d quark 
(p(k 2 ,a 2 )\d a d 1 j3 \p(ki,ai)) = 1 u p (k 2 . a 2 ) u£(&i, a x ) 



(B6) 



while for the neutron we have 

(n(k 2 ,a 2 )\u a M^|n(A;i,cri)} = 1 u^(k 2l a 2 ) u^(h, cti) 



(B7) 



(n(k 2 ,a 2 )\d a d^\n(ki,ai)) = 2 u"(fe 2 , <r 2 ) «£(&i, 01) . 

(B8) 

Using relations like [A,BC]+ = [A,B] + C - B[A,C]- 
analog equations for the four-quark condensates can be 
obtained. Two illustrative examples are given for the 
flavor- unmixed four-quark condensate for the proton 



d 3 r lC ikiri J d 3 r 2 e~ ik2r2 
</£ 2 (r 2 ,0) 



u a u l u 7 u 3 s ,ip p (n,0) 



10} 



= ~ (uu)q ^3( 70 )5 Ql (7o)q 2 7 <W 

— (7o)«ai (7o)a 2 a ^ Pi + 3 (7o)/3qi (7o)a 2 a <5 7 <5 

-(7o)/3ai (7o) (B9) 

where the normalization (|13f) and (^ttj^o — 
j2<5 1 '' S a p (uu)o [l| has been used. For four-quark 
condensates with Gell-Mann matrices involved we find 
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^p 2 (r 2 ,0) 



\ a u 3 u 7 u\,ip p (ri,0) 



d 3 n e ikl ri J d 3 r 2 e" lk2 1-2 
= |A p | 2 e abc e a ' bV ( U a ' T C 75 rf b ') (u aT C 75 d b ) ^a 11 ^-!^^ 

X (7o)<5 ai (7o)q 2 7 <W S Cl S C k 5 lj - (7o)<5q! (7o)a 2 a <W <5 Cl <5 1C S ik 



|0> (A a ) IJ (A°) 



+ (7o)/9ati (7o)a 2 a &7<5 <5 J C (5 X ° <$ k 1 - (7o)/3ai (7o)a 27 $aS S i ° 5" V 



(7o)i5ai (7o)a 2 a <W + (7o)/3ai (7o)a 2 7 <5a<5 



0- 



By using the same technique the general result for the 
six-quark condensates is obtained as 



J 



(BIO) 



-7 3 v e -'k 2 r 2 



C 2 ( r 2>o) 



"a "/3 S ^ *C ' % ( ri '°) 



\ Ap \ 2 e abc e a ' bV ( u a ' T C 75 d b ') (V T C l5 { (fij, B$ 4 C } ( 7o ) QlC ( 7 o)a 2 a S^' 1 



-{u a tip u\ u™)o (7o) ai C (7o)a 27 <5 cn <5 c k + (w^ u!^ 4)o (7o) Ql c (7o)a 2C £ cn <5' 



cncc rn 



+ (u 3 w") (7o)ai« (7o)a 2 a<5 Cl <5 C ' + (u a U e m"} (7o)ai« (7o)a 27 £ Cl # 



(U Q u")o (7o)ai<5 (7o)a 2 e ( 5 Cl ^ C m + (u~ U £ " m"} (7o)ai/3 (7o)a 2 a (5°^° ' 



cl rc'm , /-?- k -1 -?- m „-n\ 



(u a u\ ii™ v%) (70)01/9 (7o)q 27 <5 cj <5 C k - (u a " 7 4 w")o (7o)ai/9 (7o)a 2£ S C1 8' 



(Bll) 



Note that for applying the normalizations (|13f) and 114|) 
one needs a term 5 c c , which naturally arises when eval- 



uating a specific matrix element under consideration. 
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